Heat transport of electron-doped Cobaltates 
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Within the t-J model, the heat transport of electron-doped cobaltates is studied based on the 
fermion-spin theory. It is shown that the temperature dependent thermal conductivity is charac- 
terized by the low temperature peak located at a finite temperature. The thermal conductivity 
increases monotonously with increasing temperature at low temperatures T < O.IJ, and then de- 
creases with increasing temperature for higher temperatures T > O.IJ, in qualitative agreement 
with experimental result observed from Na^CoC^ . 
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The sodium cobalt oxide Na^OoO^ has become of con- 
siderable interest in the past few years because of its 
unconventional physical properties 1 . The structure of 
Na^CoO^ consists of triangular C0O2 sheets separated by 
layers of Na + ions 1,2 . This structure is similar to cuprates 
in the sense that they also have a layered structure of the 
square lattice of the Cu02 plane separated by insulating 
layers 3 . In the half-filling, the cobaltate is a Mott insu- 
lator, then the system becomes a strange metal by the 
electron doping 1 ' 2 . Moreover, superconductivity in the 
hydrated system Na^CoG^ • I.3H2O has been observed 
in a narrow range of the electron doping concentration 4 , 
around the optimal doping x ~ 0.3. Although the ferro- 
magnetic correlation is present in Na^CoC^ for the large 
electron doping concentration (x ~ 0.7) 5 , the antiferro- 
magnetic (AF) short-range spin correlation in Na^OoC^ 
and Na^CoC^ • yB^O in the low electron doping concen- 
tration (x rj 0.35) has been observed from the nuclear 
quadrupolar resonance and thcrmopower as well as other 
experimental measurements 6 . Therefore the unexpected 
finding of superconductivity in electron-doped cobaltate 
has raised the hope that it may help solve the unusual 
physics in doped cuprates. 

The heat transport is one of the basic transport prop- 
erties that provides a wealth of useful information on car- 
riers and phonons as well as their scattering processes 7 . 
In the conventional metals, the thermal conductivity con- 
tains both contributions from carriers and phonons 7 . The 
phonon contribution to the thermal conductivity is al- 
ways present in the conventional metals, while the mag- 
nitude of the carrier contribution depends on the type 
of materials because it is directly proportional to the 
free carrier density. However, the phonon contribution 
to the thermal conductivity is strongly suppressed for 
doped cuprates on a square lattice 8 . Moreover, an un- 
usual contribution to the thermal conductivity of doped 
two- leg ladder materials has been observed 9 , and this un- 
usual contribution may be due to an energy transport via 
magnetic excitations 9 . Recently, the heat transport of 
the electron-doped cobaltate Na^OoC^ has been studied 
experimentally 10 , where the behavior of the temperature 
dependent thermal conductivity is similar to the case of 



the doped two-leg ladder materials, and can not be ex- 
plained within the conventional models of phonon heat 
transport based on phonon-defect scattering or conven- 
tional phonon-electron scattering. Within the charge- 
spin separation fermion-spin theory 11 , the heat trans- 
ports of doped cuprates on the square lattice and doped 
two-leg ladder materials have been studied 12 ' 13 based on 
the t-J model, and the results are in agreement with 
experiments 8 ' 9 . Since the strong electron correlation is 
common for doped cuprates, two- leg ladder materials, 
and cobaltates, then these systems may have the same 
scattering mechanism that leads to the unusual heat 
transport. In this paper, we apply the charge-spin sep- 
aration fermion-spin approach to study the heat trans- 
port of electron-doped cobaltates. Our results show that 
the thermal conductivity of the electron-doped cobal- 
tate Na^CdC^ is characterized by the low temperature 
peak located at a finite temperature. Our results also 
show that although both dressed charge carriers and spin 
are responsible for the heat transport of electron-doped 
cobaltates, the contribution from the dressed spins are 
dominant. 

As in the doped cuprates, the two-dimensional C0O2 
plane dominates the unconventional physics of electron- 
doped cobaltates. It has been argued that the essential 
physics of the doped C0O2 plane can be described by the 
t-J model on a triangular lattice 14 , 
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where summation is over all site i, and for each i over 
its nearest- neighbor fj, C\ a (Ci a ) is the electron cre- 
ation (annihilation) operator, Si = CjaCi/2 is the spin 
operator with a = {(y x ,a y ,a z ) as the Pauli matrices, 
fi is the chemical potential, and the projection oper- 
ator P removes zero occupancy, i.e., C\ a Ci a > 1. 
In order to use the charge-spin separation fermion-spin 
transformation 11 , the t-J model (1) can be rewritten in 
terms of a particle- hole transformation d a — > f\_ a as 15 , 
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H = t'Eflfi+H* +nT,flfto + JJ2 Si- S i+ „ (2) = iheJ^f^iv ff)[S} S^ + ,,S Z _, 



then the original local constraint ^2 a Cj^Cic > 1 is trans- 
ferred as ^ufiafia < 1 to remove double occupancy, 
where f} a (fa) is the hole creation (annihilation) op- 
erator, and Sj = f\ofij1 is the spin operator in the 
hole representation. In this case, the hole operators 
can be expressed as fa — a\^S^ and fa = a^Sf in 
the charge-spin separation fermion-spin representation 11 , 
where the spinful fermion operator a ia = e~ l ® i ° en de- 
scribes the charge degree of freedom together with some 
effects of the spin configuration rearrangements due to 
the presence of the doped charge carrier itself (dressed 
charge carrier), while the spin operator Si describes the 
spin degree of freedom (dressed spin), then the local 
constraint, £ CT f} a fa = S+a^a^Sr + Srai^Sf = 

aia\(SfSf + S^S+) = 1 - a\a t < 1, is satisfied in 
analytical calculations, and the double dressed charge 
carrier occupancy, a\ IJ a\_ IJ = e 4 * iCT a|aje 4 * i_<T = and 
ai a a,i- a = e~ l * io aiaie^ 1 ^^" — 0, are ruled out auto- 
matically. It has been shown 11 that these dressed charge 
carrier and spin are gauge invariant, and in this sense, 
they are real and can be interpreted as the physical ex- 
citations. Although in common sense a,i a is not a real 
spinful fermion, it behaves like a spinful fermion. In this 
charge-spin separation fermion-spin representation, the 
low-energy behavior of the t-J model (2) can be expressed 
as, 

H = tJ2(anS+al +m Sr + . + ail Sr a\ +f)l S+ + J 

if) 

- A* X] a \e a ™ + Joff ' X Sl ' Sl+f l> ( 3 ) 
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with Jeff = (1 - S) 2 J, and 5 = (a\ a a i<T ) = 
electron doping concentration. 

According to the linear response theory 16 , the thermal 
conductivity can be obtained as, 



k(w,T) = - 



1 Imn Q (w,T) 



u> 



(4) 



with Hq(lu,T) is the heat current-current correlation 
function, and is defined as, 



U Q (r-T') = -(T T j Q (T)j Q (r')), 



(5) 



where r and r' are the imaginary times, T T is the r order 
operator, while the heat current density Jq is obtained 
within the Hamiltonian (3) by using Heisenberg's equa- 
tion of motion as 16 , 

j Q = i £ Rn [Hi, Hj] = j Q f) + j ( Q s) , (6a) 
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i ^i — fi i—fi+f} 1 ] 
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where Ri is lattice site, the spin correlation function \ = 
(S^S^ + -), e = 1 — 2t(j>/J e g, and the dressed charge carrier 

particle- hole parameter <j> — (at^a^^). Within the t-J 
model, the thermal conductivity of hole-doped cuprates 
on the square lattice has been obtained 12 . Following their 
discussions, the thermal conductivity of electron-doped 
cobaltates on the triangular lattice can be obtained as, 



k(ui,T) = K f (uj,T) + n s (u),T), 



(7a) 



1 f° 
= E A hsk J ^ -^A fa (k,w'+u) 

. ,, n F (u' + u) - n F (uj') 
xA /CT (fc,w')— ^ — , (7b) 

*Mk,u/) nB V + »l- nB V\ (7c) 

where A f = -Z X t{n + Z X t lk ), A s = (ZJ cB ) 2 e(2e X + 
2C - 4 X 7fc), l 2 sk = {[^K + sm(k x /2)cos(V3k y /2)} 2 + 
3[sin(V3fc v /2)cos(fc ;E /2)] 2 }/9, lk 
[cos(fc x ) + 2cos(fc a; /2)cos(v / 3fciy/2)]/3,the spin spectral 
function A s (k,ui) = —2lmD(k,oj), the dressed charge 
carrier spectral function Af a (k,w) = — 2lmg r7 (k, u), and 
np(uj) and ng(w) are the fermion and boson distribu- 
tion functions, respectively while the full dressed charge 
carrier and spin Green's functions are evaluated as, 



ft, 1 = 
D~ 1 {k,oj) = 
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gf ] 1 (k,w)-^ f (k,wY 



D(°)- 1 (fc,cj)-E s (fc,w)' 



(8a) 
(8b) 



respectively, where the mean-field (MF) dressed charge 
carrier Green's function g^ 1 (k,uj) = ui — the MF 
spin Green's function D^ ' _1 (fc,cj) = (lu 2 — ojfy/Bk, with 
B k = \[2 X z (ej k - 1) + x(lk - e)], A = 2ZJ eff , the spin 
correlation function X z = (S?S? + ~), and the MF dressed 
charge carrier and spin excitation spectra are given by, 



6c = -Zt X j k - {j,, 

u\ = A^kf + A 2lk + A z , 



(9a) 
(9b) 



respectively, where A\ — ae\ 2 (e X z + X /2), A2 = 
-e\ 2 {a( X * + e X /2) + {aC* + (1 - a)/(4Z) - ae X /(2Z)) + 
(aC + (1 - a)/(2Z) - a X z /2)/2], A 3 = X 2 [aC z + (1 - 
a)/{AZ) - ae X /(2Z) + e 2 (aC + (1 - a)/(2Z) - a X z /2)/2], 
the 
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spin correlation functions C 

C* = (l/£ 2 )£-,<S^S, 



i+r)' 1 



while the second-order 
dressed charge carrier and spin self-energy functions are 
evaluated by the loop expansion to the second-order as, 
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Auip'Uip+pi 
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Ft\k,p,p') 



E«(fc,w) = 
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(10a) 
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F^(k,p,p') 
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Cp' + Uk+p J 

(10b) 



respectively, where F{ l \k,p,p') = n F (£ p+k )[n B (ujp>) - 

n B {up + pi)] + n B (uJp+p')[l + n B (u) p >)], F^ h '(k,p,p') = 
n F (£ p+k ) [n B (ujp>+p) - n B (ay )] + n B (u) p > ) [1 + n B (u) p > +p )] , 

F^ h \k,p,p') = n F (t P+k )[l + 

n B (up + p<)+n B (ujpi)}+n B (ujpi)n B (ujp + p,), F[ h \k,p,p') = 
nF^p+kji 1 + n B (ujp +p >) + n B (u) p >)] - [1 + ra B (ay)][l + 
n B (uj p+ p,)}, F[ s) (k,p,p') = n F (£ p+p >)[l - n F {£ p ,)} - 

n B {uJk+p)[n F {£,p>) - n F (£p +p/ )}, and F^ s \k,p,p') = 
n F (£ p+p/ )[l - n F (£ p ,)] + [1 + ra B (w fe+p )][n F (£ p /) - 
n F (^p+ P ')). In order not to violate the sum rule of the 
correlation function (S^S^) = 1/2 in the case with- 
out AF long-range order, the important decoupling pa- 
rameter a has been introduced in the MF calculation 17 , 
which can be regarded as the vertex correction, then all 
the above MF order parameters, decoupling parameter 
a, and chemical potential fi are determined by the self- 
consistent calculation 18 . 

We are now ready to discuss the heat transport of 
electron-doped cobaltates. The observable temperature 
dependence of thermal conductivity k(T) in the exper- 
iments can be obtained from Eq. (7a) as n(T) = 
lim^^o k(lu, T). We have performed the numerical calcu- 
lation for this thermal conductivity k(T), and the results 
of k(T) at the electron doping concentration x = 0.31 
(solid line), and x — 0.33 (dotted line) for t/J = 2.5 
are shown in Fig. 1. For the comparison, the experi- 
mental result 10 taken from Na3.iCo02 is also shown in 
Fig. 1 (inset). Our results show that the thermal con- 
ductivity k(T) increases monotonously with increasing 
temperatures for the low-temperature range T < 0.1J, 
and reaches at the peak position in the temperature 
T w 0.1J, then decreases with increasing temperatures 
for the higher temperature range T > 0.1J. Using a rea- 
sonable estimation value of J sa 300 to 500fc in Na^CoC^, 



the position of the peak is at T « 0.1J « 30 ~ 50k. 
These results are in qualitative agreement with the ex- 
perimental data 10 . During the above calculation, we have 
found as in doped cuprates 12 that although both dressed 
charge carriers and spin are responsible for the ther- 
mal conductivity n(T) 7 the contribution from the dressed 
spins is much larger than these from the dressed charge 
carriers, i.e., k s (T) 3> Hf(T) in electron-doped cobal- 
tates, and therefore the thermal conductivity of electron- 
doped cobaltates is mainly determined by its dressed spin 
part k s (T). In this case, the physical interpretation of 
the present results is the same as in the case of doped 
cuprates 12 , i.e., the observed unusual thermal conduc- 
tivity of electron-doped cobaltates is closely related to 
the incommensurate spin response in the systems 19 ' 20 . 
Since n a (oj,T) in Eq. (7c) is obtained in terms of the 
dressed spin Green's function D(k,u>), while the dynam- 
ical spin structure factor of doped Mott insulators on the 
triangular lattice has been obtained from the dressed spin 
Green's function as 19,20 , 



S(k, uj) = -2[1 + n B {uj)}lmD(k, uj) 

= -2[l + n B (w)]B k 1mZ s (k,u) 

[uj 2 -uj 2 - Re£ s (fc, uj)} 2 + [Im£ s (fc,w)] 2 , 



(11) 



where Im£ s (fe,o;) and ReS s (/c,w) are corresponding 
imaginary part and real part of the dressed spin self- 
energy function Y, s (k,u)) in Eq. (10b). As we have 
shown in detail in Ref. [19,20], the dynamical spin 
structure factor (11) has a well-defined resonance char- 
acter. S(k,u) exhibits a peak when the incoming neu- 
tron energy uj is equal to the renormalized spin excitation 
E 2 = lo 2 + _BfcRcS s (fc, Ek) for certain critical wave vec- 
tors kg (positions of the incommensurate peaks), then 
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FIG. 1. The thermal conductivity k(T) as a function of 
temperature at the electron doping concentration x = 0.31 
(solid line) and x = 0.33 (dotted line) with t/J= 2.5. 
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the height of these peaks is determined by the imaginary 
part of the dressed spin self-energy function ImSs^^w). 
Since the time scale of this dynamical incommensurate 
correlation is comparable to that of lattice vibrations 
as in the case of the square lattice 12 , then this dy- 
namical spin modulations dominate the heat transport 
of electron-doped cobaltates, i.e., energy transport via 
magnetic excitations dominates the thermal conductiv- 
ity. On the other hand, the dynamical spin response is 
doping dependent, this leads to the thermal conductivity 
of electron-doped cobaltates also is doping dependent. 

In summary, we have studied the heat transport of 
electron-doped cobaltates within the t-J model. Our 
results show that the thermal conductivity of electron- 
doped cobaltates is characterized by the low temperature 
peak located at a finite temperature. The thermal con- 
ductivity increases monotonously with increasing tem- 
peratures at low temperature range T < 0.1 J, and de- 
creases with increasing temperatures at higher temper- 
ature range T > 0.1J. Our results arc in qualitative 
agreement with the experimental data of Naa;Co02. 
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